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Abstract 

We find that in generic field theories the combined effect of fluctuations and interactions leads to 
a probability distribution function which describes fractional Brownian Motion (fBM) and "com- 
plex behavior" . To show this we use the Renormalization Group as a tool to improve perturbative 
calculations, and check that beyond the classical regime of the field theory (i. e., when no fluctua- 
tions are present) the non-linearities drive the probability distribution function of the system away 
from classical Brownian Motion and into a regime which to the lowest order is that of fBM. Our 
results can be applied to systems away from equilibrium and to dynamical critical phenomena. We 
illustrate our results with two selected examples: a particle in a heat bath, and the KPZ equation. 
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Complex behavior is ubiquitous QTJ . From fluids to ecosystems to chemistry, we are famil- 
iar with phenomena often associated with non-gaussian Probability Distribution Functions 
(PDFs); phenomena described by PDFs with "long tails" or "stretched exponential" behav- 
iors 0. The PDFs associated with these phenomena have the property that, typically, |3| 
"improbable (very bad) events are much more likely than with a Gaussian" . Many of these 
phenomena are more or less loosely associated with, e. g., "complexity phase transitions" 
0] and "fractal behavior" ||; furthermore, often, "what is seen depends on the size of the 
observer" || and is accompanied by the property that "one law leads to many behaviors" 
H]. They occur in systems that are extended in space and in time, with many interacting 
components and where both, "the random and the regular" ||, are at work. 

Here we show how this "complex behavior" can be understood as a natural consequence 
of space-time evolution, fluctuations and interactions in a many body system observed 
at different scales and/or with varying initial conditions on the parameters describing the 
system. We do this for a generic many body system: a (quantum, thermal or stochastic) field 
theory where space-time evolution and interactions represent the "regular" and fluctuations 
represent the "random". 

The logic we use is the following: as is well known [J], in any field theory fluctuations 
and interactions (a) modify the n-point correlation functions of the field at different space- 
time points, and (b), lead to divergences in the original set of parameters such as couplings 
and diffusion constants, masses, etc. defining the system. From (a) it follows that the 
characteristic function of the field theory is modified from its originally normal and gaussian 
character, and therefore the PDF itself is also modified; furthermore, (b) leads to a scale- 
dependence of the parameters and correlation functions which, asymptotically, are power 
laws || with exponents calculable within the framework of the Renormalization Group (RG) 
i, 0- 

Specifically, we address here the problem of how interactions and fluctuations manifest 
themselves at asymptotic scales on the system's probability distribution and discover that 
(i) the specific asymptotic properties of the PDF depend on the basins of attraction of 
various RG fixed points and (ii) that, asymptotically, the field theory displays the type of 
behavior we have called "complex" above. 

We begin by introducing the random field <f)(x,t x ). We coarse-grain (or filter) this field 
by means of a "window-function" Wr which averages out the small scale features in 0; here 
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we mean small in comparison to a window length scale R. Without loss of generality, we 
take the window function to be translationally invariant and define 



<t>n{v) = y J dx W R (y - x) <j>(x), (1) 

where it is understood that x = (x,t x ) where x is a ^-dimensional vector, and V = flT 
is the spatial volume Q times the time interval T over which the field is filtered. Typical 
window functions are the "top-hat" window Wr(x) = Q(R — \x\) and the Gaussian window 
Wr(x) = exp(— \x\ 2 /2R 2 ). In ([[]) the coarse-grained field <pR(y) is the result of averaging 
over a space-time region or "cell" of linear extent R, "centered" at space-time point y, over 
the time interval T. The small scale features of the field are blurred within the space-time 
region selected by the window: the resolution is degraded. We also note that ([I]) is a linear 
transformation relating the microscopic (0) to the coarse-grained field (4>r) that need not 
be invertible. 



We can calculate |Tl|, the probability that the coarse-grained field 4>R(y) takes the 



value ip within the window: 

p(<p; R, y) = (5((f> R (y) ~ V>))p = J W] &(Mv) ~ <f) P \<P\ , (2) 

where P[<f>] is the PDF for the microscopic field configurations and the integral is a path 
integral over all field configurations. We point out that J dipp(ip) = f[D(f)]P[<fi\ = 1, so 
that p is normalized to unity if and only if the PDF is. We see moreover that p is a 
probability density and scales dimensionally as the inverse field y? -1 . We will make use of 
this fact below. The above probability p can be easily related to Z[J], the characteristic 
functional with source J (more commonly known as the generating functional) for the n— 
point functions. In fact, by using the Fourier integral representation of the delta distribution 
in Eq. (0) we have 

p(<p;R,y) = jdij [V<j>] e W«(w)-v]p[0] 
= J d£e~ iiv J[V(j)} e 4 1 dxW R^-y)^) p[<j)] 

/oo C 
dte-^Z[J R {x-y)=i±-WR{x-y)\ , (3) 
-oo V 
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where Z in the last line is the function of £ that results from evaluating the generating 
functional for a window-source function. For a generic gaussian probability functional [|7|], 
and |Ll 



P[0] = J^ exp !--J cfe J dz(j)(x)G~ 1 (x,z)(j)(z)\ , (4) 

the path integral in Eq.(^) is gaussian (TV is a normalization factor) and one immediately 
obtains 



where the compact notation [Wjj ■ G ■ Wr] = J dx J dz Wr(x — y)G(x, z)Wr(z — y) stands for 
the filtered or coarse-grained two-point correlation function associated with the microscopic 
field. Note that p(<p) in (||) is normalized to unity. We will calculate the asymptotic scaling 
form of [Wfl-G-Wfl] below. Here G(x, z) is the "variance" (or "dispersion") of the microscopic 
field between space-time points x and z. More generally ||13| , G(x,z) is the connected 2- 
point correlation function (or propagator), for the field 0, (cj)(x)4>(z)) = G(x, z). 

In the presence of interactions and fluctuations the characteristic functional Z[J] is mod- 
ified. It can be obtained as an expansion in terms of the n-point correlation functions for 
the field 4>(x, t), which themselves are corrected due to interactions and fluctuations. In fact, 



as is well known (T^j, the n-point functions obey Renormalization Group equations (which 
follow |0| , |TU| from the fact that removal of divergences in the model introduces an arbitrary 
scale fx) describing how n-point functions change as the parameters- the couplings {gj} and 
mass m- in the model, or the scale at which the system is observed, are modified. The solu- 
tion to the RG equations are the so-called "improved" n-point functions. For the connected 
n-point correlation function in momentum space, the RG equation in a mass-independent 
subtraction procedure is 



Gfliq; -.) = . (6) 



Here, 7^ is the anomalous dimension of 4>(x, t), given by 7^ = /i(<91n Z$jd\i) where Z^ is the 
wave function renormalization constant of 0, 5(g) = — /i(<91n Z m /dfi) where Z m is the mass 
renormalization constant, and q represents the momentum and frequency variables. 
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A RG equation for the improved probability density can be derived directly from Z as 
follows. Assuming a renormalizable field theory, the relation between the bare and renormal- 
ized generating functional is Z[J , {<?0i}> m o, A] = Z[J, {<?«}, m, (/]. The functional written in 
terms of the bare parameters and bare source (all having the zero-subscript) and cut-off A 
does not know about the arbitrary finite scale parameter \x , so that /j ^M = o. The 
chain-rule immediately implies that 

+ £ ^ + 5 ^ m -L + y / dx J{x) ike) ) m {9 * h m ' ^ ] = 0; (7) 

where the coefficient functions 

^ = m, (8) 

din 

H~q- = 8(gi)m, (9) 

describe the scale dependence for the couplings the mass m (if there is any) and the 
source function J. We emphasize that the renormalization of the source is equivalent to 
wavefunction renormalization, an important fact that is emphasized by Brown [|16|]. The 
source function acts as just another bare "parameter" of the theory, and it gets renormal- 
ized along with the other couplings. This RG equation for Z[J\ holds for arbitrary source 
functions. Due to (H) we now substitute J(x) — > yW R (x), into (|7]) and then by means of 
the identity (proof: expand Z[J] out in a functional Taylor series in powers of J), 



J dxJ(x)—Z[J] = J dxW R (x)-—-rZ[i{iW R /V} = Z-Z^W R /V}, (11) 



we obtain 



T(x) 1 J J HK J 6W R (x) 

> 1 ^^ +5{9i)m 'L Jrl i i U 



+ E Pi-^ + 6 ^)m— + f £-)z[itW R /V, { 9i }, m, /x] = 0. (12) 
Following (HD we Fourier transform this to arrive at the RG equation for pi mp : 

fair. + E Pjjrr + 6 ^ m jL ~ "T jr-^Pimpiw {</*}, m, //; R) = 0, (13) 



which follows after an integration by parts to eliminate ^-derivatives in favor of yj-derivatives. 
Just as does equation @, this equation expresses the independence of the physics on the 
choice of scale \i ~ at which we defined the values of the coupling constants gi of the 



model. The coefficient functions f3j, which are calculable using perturbation theory, describe 
the scale dependence for each of the couplings according to (§) . The scale fi is known as the 
"sliding" scale, and represents the scale at which the system is sampled. A hurricane looks 
very different depending on whether it is seen by a fly trapped inside it or by an astronaut 
from outer space; in one case fj,pi y ~ where £pi y is the typical length scale for a fly 

and correspondingly for the astronaut, where now £ Astronaut is the size of the region of the 
Earth observed by the astronaut. y?(/i) satisfies the differential equation /x ^^ = —^^(fj,) 
(compare this to fllPf) ) and represents the value of the coarse-grained field at scale fi. The 
explicit relation between /i and R is discussed in more detail below. 



The solutions to equations (0) (for n = 2) and (O) are respectively given by 



Gi mp (q; ...; gj(ji), m(//); //) = 
e L 1 ' t > {u)dlnu . G Imp (q ; ...; ^(/io), m(/i ); /i ) , (14) 

and 

1 f M 7 (u) d\nu 

P/ mp (<^(/i), {#«(//)}, m(/i),//;i2) = e 2j fo 70 " nU Pi mp (v(Vo), {9i(jk>)}, m(fj, ), fi ] R). (15) 

These solutions are to be interpreted as follows: for example, to obtain the improved 
|T7| 2-point correlation function Gj mp , one needs to write down its explicit form at some 
scale fj,Q where it is known, and the values of the various couplings gi and mass m must be 
substituted by their "running" or "effective" value gi(fJ,) and m(/i) which are the solutions 
to the RGE equations (HH). 

As one approaches a fixed point of Eq. (§), the couplings on which the anomalous 
dimension 7^ depends go to constant values g*, and the anomalous dimension reaches a 
constant value; similarly, as a result of equation (|6]), the two-point correlation function 
Gj mp (q; ...; g(fi); fj) goes (after Fourier transforming) into a function f{u) of argument, u = 
t/r z (where z is the value at the fixed point of a different anomalous dimension, the so 
called dynamical exponent) times r 2x , where \ is related to 7^ at the fixed point. Here, 
r = \x — z\, and t = \t x — t z \. In the neighborhood of any RG fixed point, it is easy to 
derive the scaling form of the Green functions of the field theory in terms of dynamic critical 
exponents \i z as follows. Under independent rescaling of coordinates and the time x = sx', 
t = s z t' (note: this kind of independent space and time scaling is needed when treating 
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non-relativistic theories such as arise for example in diffusion and growth processes) the 
field scales as <fi(x,t) = s x (j)(x' ,t'), so it follows that 

G(x,t) = (0(f,t)0(O,O)> 

= s 2x {(j){s' l x, s-^)0(O,O)), 

= s 2x G{s~ l x , s~ z t) 

= r*f(^), (16) 

where the last line follows from choosing s ~ r = \x\. The asymptotic behavior of f(u) 
is given by Jl8 



lim f (u) — > const. , and lim f(u) — > u 2x ^ z . (17) 

To aid our understanding of the relation between the coarse-graining scale R and the 
sliding (momentum) scale fi, we can use the above general solution ([15D plus simple dimen- 
sional analysis. Let d v denote the canonical dimension of the field: [</?] = then p(<p) 
must have dimension [p] = [y? -1 ] = yT dv expressed in units of the sliding scale [i. The mass 
dimension [m] = ft, and assuming dimensionless couplings, then [gj] = fi° = 1. Then from 
(|15|) it is easy to prove that 

^(^oH^o)},^^ (18) 



where F is a dimensionless function of the dimensionless arguments as wriiten here. Cast in 
this form, we can investigate the infrared limit of the probability density by taking /i — > 0. 
The connection to the Wilsonian RG approach and the window scale is the following. In 
the Wilsonian RG, the degrees of freedom are coarse-grained either in real space or in 
momentum space, the latter typically proving to be the more technically convenient choice. 
This is performed over a finite- width momentum shell corresponding to A/s < \k\ < A; 
this corresponds in fact to a "top-hat" window in momentum space. The UV cutoff A is 
contracted down to A/s, where s > 1. The infrared limit obtains by taking s — ► oo. The 



connection with the field theory sliding scale is fi = A/s []19 |. Therefore, when coarse- 
graining in real space, the contracted cutoff corresponds to an increasing length scale: R = 
s/A. So, the IR limit fj, — > corresponds to R — > oo. The dimensionless product fiR = 1, 
and we can now replace /i = 1/R everywhere in (p!8|). In the neighborhood of an IR fixed 
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point, reached by taking R — >• oo, we therefore obtain the asymptotic scaling form: 
lim pjmpitpQkddQiMh rn{n ) , fj, ; R) = B?* + b*W) 

x F /^ ( iL ) ^-h^) j !^(iL)*to-M ) rf) ) (19) 

v no* fMy Ho yUo J/ 

where /io = l/-Ro is some reference scale and we have used the solutions of (|8]||) at the fixed 
point g* in arriving at this final form. At the fixed point, the exponential prefactor in flTS| ) 
scales as (////io) -57 ^ 9 ^ ■ 

Note that ( plf) demonstrates in general that the large scale asymptotic form of the coarse- 
grained probability density goes as a non-trivial power of the window size R times a certain 
dimensionless function. The asymptotic behavior is controlled by the canonical dimension 
d^, of the field, its anomalous dimension at the IR fixed point 7^(5'*) (as well as by S(g*) 
in a massive theory). Specification of pimp at some reference scale Rq = l///o yields the 
explicit mathematical form of F. The derivation of the RG equation for pim P in flUf ), and 
the large-distance scaling behavior of its general solution in fll9| ) are the key results of this 
Letter. 

To illustrate the use of the above general results, we apply them to two simple examples: 
(i) a free particle in the presence of a heat bath, and (ii) a system described by the Kardar- 
Parisi-Zhang (KPZ) equation with colored noise |2(| . 

In the case of a free particle in a heat bath there is no self-interaction, but the statistics of 
the bath turns the problem into a classical Brownian motion problem; for the KPZ system 
there are interactions among the particles making up the system (described by the KPZ 
field) and also interactions with the bath (which could be an external environment or the 
effective result of a "microscopic" dynamics) represented by a noise term that drives the 
time derivative of the KPZ field. For each of these examples we will obtain the form of 



equation (|T8[) which corresponds to the scale-dependent form of the probability pi mp ((p; R). 
To lowest order, the probability density is gaussian, thus the procedure consists in computing 
the improved form of G(x, z; t x , t z ) appearing in Eq. and given by Eq. (|14|). We indicate 
briefly how to RG-improve non-gaussian probabilities below. 

Example (i). The equation of motion for a point particle (a particle may be regarded as 
a zero-dimensional field) in a medium with viscous drag and subject to a random force is 
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(20) 



Here P of Eq. (|) is |4 



1 /•* , /cF x 2 



PHcexp^-^rf^-j > , (21) 

where f(t) denotes the coordinate of the particle at time t and v (t) is the particle's velocity. 
The object G(t, t') is 

Gij(t, = (r^r^t')} « t 6^6(1; - t') . (22) 

There are no corrections due to fluctuations or interactions, and the probability of Eq. 
(|) is simply given by inserting Eq. (^) into (H) (after using a temporal "top-hat" window 
W T ) 

K r ) a ^' ex p(-^) > ( 23 ) 

(where d is the number of components of the vector r) which of course is the probability 
distribution that a particle executing standard Brownian motion be at position r = \r\ at 
time T. 

Example (ii). The KPZ equation (|24] ) is a non-linear Langevin equation for a field. 
Contrast this to (|20|) which is a linear Langevin equation for a point-particle. In a system 
described by the KPZ equation 

+ ^ = ^V 2 0+^A(V0) 2 + r ? (x,t) (24) 
with colored noise |2l[] r/(x, t), there are corrections due to both fluctuations and interactions. 



Setting A = yields the linear Edwards- Wilkinson (EW) model having a unique IR fixed 
point PI, for which the IR critical exponents are known exactly for all space dimensions: 
(z,x) — (2,(2 — d)/2), in the case of white uncorrelated noise [p2]1 - The EW model is 
the free-field limit of the KPZ equation and generalizes the concept of random walk of a 
classical particle to the level of free fields. For non-zero A two new fixed points arise and 
their corresponding critical exponents x an d z have values differing from the EW fixed 
point exponents. Briefly, in d = 3 space dimensions PI is a saddle point with exponents 
(z,x) = (2, — |), P2 is infrared unstable with (z,x) = an d -P3 is infrared stable 

with (z, x) — (Ij !)■ For each of these fixed points one has to consider two possibilities (cf. 
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Eq. above), depending on whether \t — t'\ ^> \f — r*\ z or \t — t'\ <C \r — r*\ z since, as 

mentioned above, these limits lead to different asymptotic behaviors for the scaling function 
/ in (|T7|). The corrected form of the PDF that the effective field has value (p is then obtained 
by inserting the "improved" form of G, Eq. flTBp, into Eq.(^); coarse-graining it with the 
window function W R and using (JTSj) leads to 



Pi mp (<f] R) 



e -5^ o T0( u ) dlnu 



V 2 ( ^ 2 

■ exp 



\2-K[W R -G Imp -W R ] l \ ^[W R -G Imp -W R }J- 
To proceed further we need the scaling form of the coarse-grained improved two-point 
function. To this end, we take a window function of the simple form W Ri t{x) = Q(R — 
\x\)Q(T — t x ) and without loss of generality, take the window center at the origin. (The 
scaling cannot and does not depend on where the window is located.) Then we find that 

^[W r ,t ■ G Imp ■ W R , T ] ~ R 2x , (25) 



for \t\ -C \r\ z and 



^[W r ,t ■ G Imp ■ W r ,t] ~ T 2 ^ z , (26) 



for \t\ 3> \r\ z , respectively. 

For each of the fixed points the asymptotic limits are 

^ Plm MR,T) ~ r-x/'+h.OrV-e^ ^ (27) 

and 

|4 gg,.P/mp(w ^ T) ~ J2-x+^&T) exp (28) 

The roughness exponent and the anomalous dimension are related through the exponent 
identity \ = —d^ — ^^(g*) which follows from using ([L4]) plus dimensional analysis to arrive 
at a scaling form for Gj mv in complete analogy to what we worked out above for pimp in 
(|18[) and in (|19|). Comparing the scaling form so obtained with fll6|) immediately yields this 
identity. 

We see that in the KPZ problem the asymptotic probability distributions associated with 



the fixed points P2, and P3 are the ones for fractal Brownian motion [23], unlike in the 
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free field case, A = 0, or EW model, which describes Brownian motion for non-interacting 
fields. 

In comparing KPZ with EW results, we have fractal Brownian motion in time, Eq (p7|). 
and in space Eq. because the exponent combinations 1\j z and 2\ appearing within 

the exponential function are different from the EW values at the fixed points of the RGE 
P2 and P3. 

Moreover, a nonzero wavefunction renormalization 7</,(g*) 7^ can modify the exponents 
of the power-law prefactors in ([27]) and (^8|) away from their naive canonical values. For 
EW, the anomalous dimension is identically zero. The reason for the appearance of fractal 
Brownian motion in Example (ii) is now obvious: the combined effect of fluctuations and 
interactions drives \ and z away from their free field theory (EW) values and the probability 
distribution is correspondingly shifted from regular Brownian motion to fractal Brownian 
motion. 

RG-improvement and the asymptotic scaling of non-gaussian probabilities can be worked 
out using the techniques of this paper. For any interacting field theory, one can expand the 
exact characteristic functional about the gaussian limit and proceed to derive the associated 
RG equation for pi mp by substituting the identity 

Z[J] = exp (^S 'inter action[-^~j]j ^Gaussian [J] \ j^i£\y R 5 (29) 

directly into the definition of pi mp in (0)- Here, Si nte raction represents the non-Gaussian part 
of the action. In our examples, we have calculated the lowest order correction, which is the 
RG improved Gaussian. Non-Gaussian terms will appear at higher order in the couplings 
and these will take the form of a polynomial in (p times a gaussian. The polynomial will 
depend on higher n-point Green functions (e.g., for n > 2), each of which obeys Eq(|6]). 

In summary we see that the combined action in space-time evolution of interactions 
and fluctuations leads to probability distributions which at some scale may be gaussian. 
However, given a particular dynamics, as the spatial extent of the system or the time window 
over which the system is observed, is changed, the two-point correlation function acquires 
an anomalous dimension and the effective PDF becomes the one for fractional Brownian 
Motion (fBM). This fBM is commonly associated with complex behavior, such as known 
to occur in the financial markets, ecology or in river systems |TJ], and its character depends 
on the initial values assigned to the relevant couplings at some scale. We have developed 
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the general framework and carried out the above analysis for two complementary problems. 
For the free field, as expected, we recover Brownian motion; this is useful as a very simple 
test of the validity of the application of our framework; in the much more complex case of 
KPZ we find a wealth of behaviors, with persistence (related to superdif fusive processes) 
antipersistence (related to suhdif fusive processes) and regular Brownian Motion P^fl . In 
fact, our analysis explicitly shows that the nature of the fBM displayed by the field system is 
related to the choice of initial conditions for the couplings, since the values of the roughness 
and dynamical exponents x an d z depend on which basin of attraction of the fixed points 
the initial couplings are chosen. Thus, the notion of complexity is scale dependent as well as 
dependent on the initial conditions; furthermore, since a given dynamics may have a variety 
of fixed points, we see explicitly that a particular system may display various complex 
behaviors which reflect the presence of different environments. 

Finally, we note that the methods presented here can be extended and applied to any 
system which can be described by a field theory, such as fluids, materials, aggregates, and 
so forth. 

Acknowledgments 

The authors thank M. Visser for discussions. J. P-M thanks M. Gell-Mann and G. West 
for many stimulating conversations on the emergence of power laws in science. 



[1] See, e.g., special issue of Science on "Complex Systems", Science 284, 2 April 1999. 
[2] H. Scher, M. Shlesinger and J. T. Bendler, Physics Today 44 (1991) 26. 
[3] N. Goldenfeld and L.P. Kadanoff, Science 284 (1999) 87. 

[4] S A. Kauffman, The Origins of Order: Self-Organization and Selection in Evolution, Oxford 

University Press, New York, 1993. 
[5] B.B. Mandelbrot, The Fractal Geometry of Nature, Freeman, San Francisco, 1982. 
[6] M. Gell-Mann, The Random and the Regular, Santa Fe Institute's 1999 Ulam's Lecture Series, 

Santa Fe, NM, November 1999. 
[7] S. Weinberg, The Quantum Theory of Fields I & II, Cambridge University Press, Cambridge, 



12 



[8 

[9 
[10 



[11 
[12 

[13 
[14 

[15 

[16 

[i7; 



England, 1996. 

K. Wilson and M.E. Fisher, Phys. Rev. Lett. 28 (1972) 248. 
M. Gell-Mann and F.E. Low, Phys. Rev. 95 (1954) 1300. 

N. Goldenfeld, Lectures on Phase Transitions and the Renormalization Group, Addison- 
Wesley, New York, 1992. 

C. Itzykson and J. B. Zuber, Quantum Field Theory, McGraw-Hill, New York 1979. 

R. J. Rivers, Path Integral Methods in Quantum Field Theory, Cambridge University Press, 
Cambridge 1987. 

D. Hochberg, C. Molina-Paris, J. Perez-Mercader and M. Visser, Phys. Rev. E60 (1999) 6343. 
J. Zinn-Justin, Quantum Field Theory and Critical Phenomena, Oxford University Press, 
Oxford, 1996, 3rd edition. 

D. Amit, Field Theory, the Renormalization Group and Critical Phenomena, World Scientific, 
Singapore, 1984, 2nd edition. 

L.S. Brown, Quantum Field Theory, Cambridge University Press, Cambridge, 1992. 



The term "improved" means ||, [0] and [10], that the effects of interactions and fluctuations 
have been taken into account. This is so on account of the fact that the coefficients in equations 
@ and ([H]) are already calculated at each order of perturbation theory. The objects pimp 

(2) 

and G Ir ^ p contain much more information. 
[18] E. Frey and U. C. Tauber, Phys. Rev. E50 (1994) 1024. 

[19] M. Le Bellac, Quantum and Statistical Field Theory, Oxford University Press, Oxford, 1991; 
Chap 7. 

[20] E. Medina, T. Hwa, M. Kardar and Y.-C. Zhang, Phys. Rev. A 39 (1989) 3053. 

[21] The noise is defined to be gaussian, with zero mean and a two-point correlation function given 
by (n(x,t)rj(x' ,t')) = 2D(x, x';t,t'). The local function D(x,x';t,t') is the double Fourier 
transform of D(k,u>) oc k~ 2p oj~ 2S . The existence of power law behavior for the correlation 
functions of eft is not due to the choice of noise correlations given above; it is a generic 
prediction of the Renormalization Group due to the fact that the RGEs are quasi-linear first 
order partial differential equations. 

[22] A.-L. Barabasi and H.E. Stanley, Fractal Concepts in Surface Growth, Cambridge University 
Press, Cambridge, 1995. 

[23] K. Falconer, The Geometry of Fractal Sets, Cambridge University Press, Cambridge, 1986. 

13 



M. Shlesinger, G. Zaslavsky and U. Frisch, Levy Flights and Related Topics in Physics, Lecture 
Notes in Physics, 1995, Springer-Verlag. 



14 



